Hyperfine spectra of trapped Bosons in optical lattices 
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We calculate the interaction induced inhomogeneous broadening of spectral lines in a trapped 
Bose gas as a function of the depth of a three-dimensional cubic optical lattice. As observed in 
recent experiments, we find that the terraced "wedding-cake" structure of Mott plateaus splits the 
spectrum into a series of discrete peaks. The spectra are extremely sensitive to density corrugations 
and trap anharmonicities. For example, even when the majority of the cloud is superfluid the 
spectrum displays discrete peaks. 

PACS numbers: 03.75.Lm, 32.30.Bv, 03.75.Hh 



I. INTRODUCTION 

The study of quantum degenerate atoms confined to 
periodic potentials forms an important subfield of mod- 
ern atomic physics. Research in this area is driven by 
its connection to condensed matter physics and quantum 
information processing [3, [3]. A rich set of probes, includ- 
ing optical spectroscopy, noise spectroscopy, interference, 
and density profile measurements 0, 0, H, @, @ > have 
been used to characterize these systems, with a focus on 
understanding the interaction driven superfluid-insulator 
transition. Here we analyze in detail what information 
one gains from inhomogeneous pressure shifts of spectral 
lines in a gas of bosons confined to an optical lattice. 

Atomic interactions lead to pressure and density de- 
pendent shifts of atomic lines. These "clock shifts" limit 
the accuracy of atomic clocks. In an inhomogeneous sys- 
tem they are spatially dependent, yielding a broadened 
spectrum whose structure reveals details about the local 
atomic correlations. In many situations the clock shift is 
proportional to the atomic density, and the spectral line 
directly gives a histogram of the atomic density. As an 
example of this technique, Bose-Einstein condensation in 
spin polarized atomic hydrogen was detected through the 
line shape of a two-photon ls-2s transition [9( . More re- 
cently, Campbell et al. Q utilized atomic clock shifts to 
experimentally probe bosons trapped in an optical lat- 
tice, finding evidence for Mott insulating shells. Moti- 
vated by these latter experiments, we present a theoreti- 
cal analysis of the lineshapes which should be found when 
bosonic atoms in a periodic potential are confined by a 
nominally harmonic potential. 

In Sec. |TT] we use a local density approximation to cal- 
culate the spectrum of a harmonically trapped gas as 
a function of the depth of an optical lattice (Fig. [1]). 
In the deep lattice limit, the spectral line splits into 
several distinct peaks, associated with the formation of 
density plateaus. Due to the sensitivity of these spec- 
tra to small density corrugations, this splitting occurs 
even when large sections of the cloud are superfluid. De- 
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FIG. 1: Theoretical spectra showing the number of 
87 Rb atoms transferred from hyperfine state \a) = 
\F — 1, Tfip = —1) to state |6) = \F = 2, = 1} when ex- 
cited by light detuned from resonance by the frequency 5v. 
The N = 9 x 10 4 atoms are confined by a harmonic potential 
with uj = (uixLUyUiz) 1 ^ 3 = 2tt x 100Hz and a three-dimensional 
periodic potential with lattice depth Vo = 5, 10, 25, 35-B rcc 
(from top to bottom). 



spite qualitative agreement with experiments, our calcu- 
lation severely underestimates the small detuning spec- 
tral weight. In Sec. IIIII we show that these discrepancies 
are consistent with trap anharmonicities. We also explore 
other sources for the discrepancy. 

Experimental Details. Since we are largely con- 
cerned with the experiment in Ref. Q, we give a 
brief review of the important experimental details. In 
these experiments, a gas of 87 Rb atoms in the \a) = 
\F = 1, tyi f = —1) hyperfine state (F the total spin and 
mp its z component) was cooled well below the con- 
densation temperature (ic| . By combining optical and 
magnetic fields, a three-dimensional periodic potential 
Vpcr = —Vo [cos(2irx / d) + cos(2iry / d) + cos(27rz/d)] was 
superimposed on a trapping potential. The spacing be- 
tween lattice sites, d = X/2 — 532nm, is half of the lattice 
lasers' wavelength. The lattice depth Vo was tuned from 

zero to 40-E r ec where E ICC = §=r (tt) * s ^ ne photon re- 
coil energy. A microwave and RF field were tuned near 
resonance for a two photon transition from the la) state 
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to an excited hyperfine state \b) = \F = 2, uif = 1) 



II. SPECTRUM OF HARMONICALLY 
TRAPPED GAS 

A. Hamiltonian and approximations 

Hamiltonian. Bosons in a sufficiently deep optical 
lattice are described by the Bose- Hubbard model [TTI ] . 
found by projecting the full Hamiltonian onto the lowest 
Bloch band. We will work with a two-internal state Bose- 
Hubbard Hamiltonian, where ai and hi annihilate bosons 
at site i in states \a) and respectively. Including an 
external trapping potential the Hamiltonian is 



<i,j> i 



2 ^i,a(^z,a 1) ~l~ ^,a^i,a 



<«,j> 



n i; b(n it b - 1) + Vj^n^i, 



^afe n i,a n i,b + H.C. 



(1) 



where n, 



def. 



The t a 's describe hopping 
rates and [/ aj g the interaction where a and /3 label 
the species (\a) or We abbreviate U a — U aa . 

We have absorbed the chemical potentials into the 
trapping potential, writing V^ a = Vj, a — /x Q where 
Vi, Q is the external potential at site i for species a. 
In terms of microscopic quantities, these parameters 
are t a = J drw* a (r) [-h 2 /(2m a )V 2 + V peI (r)] w a (r), 
Vi, a « V tla , p (Ri), and U a p = (4irh 2 a a( 3/m) J dr 
\w a (r)\ 2 \wp(r)\ 2 where m a is the mass and w a the nor- 
malized Wannier function for state a, while a a p denotes 
the a-(3 scattering length. For 87 Rb, the relevant scat- 
tering lengths are a aa = 5.32nm and a Q ;,=5.19nm [l2T |. 
The | b) atoms will be sufficiently dilute that abb will not 
enter our calculation. The competition between the ki- 
netic and interaction terms drives the Mott insulator to 
superfluid phase transition. 

In the experiments of interest, the atoms all begin in 
the | a) state, and one measures the rate at which atoms 
are transferred to the \b) state under the influence of a 
weak probe of the form i? pro be oc j bj(ije~ tuJt + H.c, 
within the rotating wave approximation, where lu is the 
frequency of the photons. To calculate this response, it 
is sufficient to understand the properties of the single- 
component Bose-Hubbard model (the terms in Eq. ([1]) 
containing only a's). 

Mean-field theory. The ground state of the single 
component Bose-Hubbard model is well approximated 
by the Gutzwiller mean-field theory (GMF) of Ref pi ]. 
This approach is exact in infinite dimensions and in the 
deep Mott insulator and superfluid limits. Sophisticated 
numerical calculations, some with a trapping potential, 



have shown that this mean field theory yields qualita- 
tively accurate phase diagrams, energies, and spatial den- 
sity profiles 0, [H, [H, [13, El. As a point of refer- 
ence, Monte-Carlo calculations predict that for unity fill- 
ing the 3D Bose-Hubbard model on a cubic lattice has 
an insulator-superfluid transition at t/U = 0.03408(2), 
while mean field theory gives t/U = 0.029. We will 
work within this approximation. As will be apparent, 
one could extend our results to include fluctuation effects 
by numerically calculating the density and compressibil- 
ity of the homogeneous system. Within the local density 
approximation discussed below these homogeneous quan- 
tities are the only theoretical input needed to determine 
the spectrum. 

The Gutzwiller mean field approximation to the 
Bose-Hubbard model can be developed either from a 
mean- field or variational standpoint. As a variational 
ansatz, GMF corresponds to taking the wavefunction 
to be a tensor product of states at each site: = 

0i (j2 n f^ K) 71 )) where \i,n) is the state with n par- 
ticles at the i'th site; the /„ are varied. In the corre- 
sponding mean-field language, fluctuations of the annihi- 
lation operators from their expectations are assumed not 
to affect neighboring sites. Then, assuming translation 
symmetry remains unbroken and letting q be the number 
of nearest neighbors, one has 



Hmf = 



-qtaUa 



n 2 



Viiii + H.c. 



(2) 



from which one can find a self-consistent ground state 
with (a) = J^n Vn + lf*fn+i- 

Local density approximation. We use a local den- 
sity (or Thomas-Fermi) approximation (LDA) to calcu- 
late the spatial dependence of thermodynamic quanti- 
ties: all physical quantities at location r are taken to 
be those of a homogeneous system at a chemical poten- 
tial ix — V(r). This is expected to be valid when the 
spatial correlation length of the homogeneous system is 
much smaller than the length scale of the trapping po- 
tential (ll. The validity of the GMF+LDA is discussed 
in Ref. 'SPJ, [H d [T3, d • 



B. Homogeneous clock shifts 

The clock shift is a density-dependent shift in the en- 
ergy splitting A for driving a transition from internal 
atomic states \a) to states \b) due to the inter-particle 
interactions. In this section we review the known results 
for the clock shift of a homogeneous system in terms of 
local correlations, and specialize to the case of atoms in 
a periodic potential. 

We will assume that t a = tb and Vi.b — V^ a + A 
where Ao is the energy splitting of the two states in vac- 
uum. These assumptions are justified in the recent ex- 
periments, where the polarizabilities of the two internal 
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states are nearly indistinguishable. In linear response, 
the average clock shift energy of the homogeneous sys- 
tem is then 



5E 2 = (U ab ~U a ) 



E 2 4 a \ a i a i 



(3) 



where the expectation is in the initial, all-|a) state [20|, 
[2lj . This expression can be rewritten in a somewhat 
more familiar form by defining the local second order 

correlation function ^2 = i^ a \ a \ a i a i^ / (^ a \ a i^ 



so 



SE, 



(Uab - U a ) g 2 (n) 



(4) 



Special cases of the clock-shift formula: dilute 
superfluid, Mott insulator, and normal fluid. For 

a dilute superfluid, the initial state is a coherent state, 
and Eq. ([3]) gives a shift proportional to the occupation 
of each site, 



SE : 



SF 



(Uab - U a )n. 



(5) 



Deep within the Mott insulating phase, the initial state 
is a number eigenstate and 



SE- 



ui 



(U ab -U a )(n-l). 



This latter formula has an intuitive explanation. In a 
Mott insulator with filling of one particle per site, the 
atoms are isolated so there is no interaction between par- 
ticles. Hence SEmi must vanish when n = 1. Fig. [5] il- 
lustrates how the clock shift energy evolves from being 
proportional to n to n — 1 by juxtaposing the contours 
of fixed 8E2 and those of fixed density. 

If one raises the temperature the system becomes a 
normal fluid, even at weak interactions. In the absence 
of interactions, the normal fluid statistical factor gi ap- 
pearing in Eq. ((4]) is g 2 — 2, [22] so that the clock shift 
energy is twice as large as in the superfluid: 



SE NF = 2(U ab -U a )n. 



(6) 



Given that there is no phase transition between the 
zero temperature Mott insulator and the normal gas, it 
is interesting that the clock shift energy changes from 
2(U a b — U a )n in the normal fluid to (U a b ~ U a )(n — 1) 
in the Mott insulator. A quantitative understanding of 
this crossover would require calculating the temperature 
dependence of the pair correlations in the strongly inter- 
acting limit. 



C. Calculation of spectrum in a trap 

To calculate the spectrum we assume that the gas can 
be treated as locally homogeneous, and we can indepen- 
dently sum the spectrum from each region in the cloud. 
As in the experiment, we imagine applying a weak probe 



at frequency co for a time r f» 100ms. Assuming that 
the finite probe duration is the primary source of broad- 
ening, second order perturbation theory implies that the 
number of atoms of atoms transferred to the \b) state will 
be 



(7) 



N b (u)(x / d 6 rn(r)S 1/T (A(r) -u), 



where S^(v) has a peak of width 7 at v = 0. We will 
model S^(v) = (l/ir)j/(v 2 + j 2 ) as a Lorentzian of width 
7. The local density n(r) and clock shift A(r) are set 
equal to that of a homogeneous system with chemical po- 
tential n(r) = jiQ — Vtrap(f)- Experimentally, the number 
of |6)-atoms is measured by monitoring the absorption of 
a laser which transfers \b) atoms into a third state. 

We calculate the integral in Eq. (0 within the 
Gutzwiller mean field approximation to the Bosc- 
Hubbard model. As illustrated in Fig. [21 both the density 
n and the clock shift A can be expressed as functions of 
the parameters fi/U and t/U. Within the local density 
approximation, t is constant throughout the trap, and fi 
varies in space, taking its maximal value /iq at the cen- 
ter of the trap. For extreme values of t/U (either large 
or small) we can analytically calculate the contours in 
Fig. 2(a) Generically, however we must rely on numeri- 



cal methods. 

Our results are shown in Fig. [1] for a harmonic trap 
Vtrap(^) = mu) 2 x x 2 /2 + muj 2 y 2 /2 + mui 2 z 2 /2. One sees 
that in the deep Mott limit, the spectrum displays sharp 
peaks, while in the deep superfluid limit, the lineshape is 
smooth. The peaks are due to the stepwise variation of 



A(/i) (illustrated in Fig. 2(a) I which lead to large regions 
of the trap where A(r) takes on discrete values. Com- 
pared with the experiments in Ref. 0] , our spectral lines 
have severely reduced small detuning spectral weight. In 
Sec. IIIII we show that trap anharmonicities can account 
for this difference. 

Note that within the local density approximation, 
the spectrum is independent of trap anisotropies, 
as long as the trap is harmonic. This generic 
feature of the LDA is seen by examining an in- 
tegral of the form I = J d 3 r /(//)=/ d 3 r f(fio ~ 
moj 2 x 2 /2— moj 2 y 2 /2— muj 2 z 2 /2). Rescaling the co- 
ordinates so that muj 2 x 2 /2 — ^qx 2 (and simi- 
larly for y and z), this integral becomes I = 



^^l/m 3 ujluj 2 uj 2 z A'K J drr 2 f(fj, a - ^i r 2 ), 



where 



V* 2 



y- 



From this analysis it is clear that apart 



from an overall scale factor, the spectral lineshape is only 
a function of the central chemical potential fig. 

Experimental Parameters. The experimental con- 
trol parameters are the optical lattice depth Vb, the num- 
ber of particles N, and the trap frequencies uj v . The nat- 
ural theoretical parameters are t, U, and ^o- To compare 
our results to experiment, we use a non-interacting band 
structure calculation to relate t and U to Vq [HI]. To re- 
late /Iq to experimental parameters we note that within 
the LDA the number of trapped atoms N is only a func- 
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tion of t/U, and fio/huj where w 3 



U} X LQ V LQ Z 



For each 



value of t/U we compute N(fi /fuli), for several values 
of no, then invert the function to get fig as a function 
of N. We then have the ability to select the value of 
/io corresponding to the number of particles used in the 
experiment. 




p 

3; 
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t/U 

(a) 





0.00 0.01 0.02 0.03 0.04 0.05 
t/U 

(b) 



FIG. 2: (a) Gutzwiller Mean Field phase diagram with con- 
stant density contours. The vertical dashed lines show the 
spatial dependance of the chemical potential for a trapped gas 
in the LDA: from left-to-right these correspond to the deep 
Mott li mit, t he "corrugated superfluid" situation appropriate 
to Fig. 4(b) and the dilute superfluid. (b) Phase diagram 
with contours of constant {71)52 = (n(n — 1)) / (n). Contours 
are spaced by 0.1, with additional lines at m±0.01, for integer 
m, to emphasize the Mott regions. 

Campbell et al. 0] do not report the number of parti- 
cles in the experiment. For Fig.[TJ we choose N = 9 x 10 4 
so that the maximum site filling for Vo — 35E lcc and 
Vo = 25E lcc is n = 5, as is observed in the experiment. 



D. Analytic Results 

1. Dilute superfluid 



where n(0) = (/io — e (fy)/U a is the central density. Sim- 
ilar expressions were obtained in Ref. [23| 

At fixed central chemical potential (equivalently, fixed 
central density) the width of the spectrum is proportional 
to U ab —U a . Given a fixed number of particles, the central 
density varies as n(0) ~ I/L/ 3 / 5 , so the width of the 
spectral line varies as Un(0) ~ U 2 ^ 5 . 



2. Deep Mott limit 

Now we analytically calculate the spectrum in the 
deep Mott insulator limit, where the density of the ho- 
mogeneous system with chemical potential /i equals the 
largest integer bounded by fi/U, denoted \n/U~\ [lHlls|. 
In the local density approximation the density jumps 
from density n — 1 to n as one moves through the lo- 
cation in space where local chemical potential is given by 
/t = U a {n — 1). Each plateau of fixed n gives a (broad- 
ened) delta-function contribution to the line shape at de- 
tuning 5u> n = 2(n — l)(U a b ~ U a ). The magnitude of the 
delta function is proportional to the number of particles 
in the plateau, leading to a spectrum 



N b (8co) oc }^ A n 8i/ T (Sui — Suj n ) 

71=1 

An = [/io - U a (n - 1)] 3/2 n 



(10) 



A n ^n = (/«o - U a (n - l)f /2 - (fjt - U a nf' * n, 

where n = \[j,o/U a ~\ is the central density. 

The deep superfluid and deep Mott insulator spectra 
are plotted in Fig. [3] using Eq.'s © and {ID]). Note the 
envelope of the spikes seen in the insulating state has 
the same shape as the superfluid spectrum. This can be 
understood from noting that in both cases the density is 
proportional to (i or \/i/U], resulting in similar coarse- 
grained A. 
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Having numerically calculated the spectra, we now spe- 
cialize to the dilute superfluid limit where the line shape 
can be calculated analytically. The clock shift energy in 
this limit is A = (U a b ~ U a ) (n), and within the local 
density approximation the site filling at position r is the 
greater of zero and 



n(r) = (mo - e(0) - F trap (r)) /U a 



(8) 



where as previously stated, Vt iap (r) — muj x x 2 /2 + 
mu> 2 /2 + muj z z 2 /2 is the trapping potential, /io is the 
central chemical potential, and e(0) is the energy of the 
k = single particle state in the periodic potential. In 
the tight binding limit, e(0) — —qt where q is the num- 
ber of nearest neighbors. Substituting this result into 
Eq. ([7]), and neglecting the broadening one finds 



N b (5uj) cx Suj^(U ab - U a )n(0) - 5uj, 



(9) 



E. Intermediate Coupling 

Finally, let's consider how the spectrum evolves as one 
increases t/U from zero. For non-zero t/U, superfluid 
shells form between Mott plateaus. These regions make 
the density continuous. Consequently, in the spectra, 
the areas of zero signal between peaks begin to fill in. 
Using our numerics, we find that the peaks remain vis- 
ible until the system is well into the superfluid regime. 
An example is shown in Fig. |4(b)[ corresponding to the 
chemical potential trajectory at t/U a — 0.018 shown in 
Fig. 2(a) Note that although the only Mott lobe crossed 
is at n = 1, six peaks are clearly visible. Clearly one 
must be cautious about using such spectra to distinguish 
superfluid and Mott insulating states. 

The source of the peaks are weak density corrugations 
which arise in the superfluid state near the Mott bound- 
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FIG. 3: Analytically calculated spectra for the harmonically 
trapped system in the deep Mott limit (solid line) , plotted as 
a function of 5u)/(U a b — U a ). Superimposed is the spectrum 
of the superfluid (dashed line) with the same parameters, but 
horizontally shifted to the left by — 8uo/2{U a b — U a ). The 
central density is n s (0) = [fio— e(Q)]/U a = 5.8. This illustrates 
that the envelope of the spectral line in the Mott insulating 
state has the same shape as the superfluid spectrum, shifted 
horizontally. 



the Mott insulator limit, a smooth distribution in the 
deep superfluid limit. Furthermore, the overall energy 
scales of our spectra are consistent with those found ex- 
perimentally. We caution however that we have used 
identical trapping frequencies in each of our spectra, 
while experimentally the harmonic confinement varies in 
an uncharacterized manner when the optical lattice in- 
tensity is changed. With this systematic variation makes 
quantitative comparison difficult. 

Despite the qualitative similarities between theory and 
experiment, serious discrepancy remains. In particular, 
the experiment finds much more spectral weight at small 
detunings than theory predicts. Here we explore possible 
sources of this discrepancy. Our primary result is that the 
discrepencies are consistent with trap anharmonicities. 

In Sec. IIII Al we give an analysis of trap anharmonic- 
ities. In the following sections we briefly discuss several 
other possible explanations of the discrepencies: non- 
equilibrium effects and nonlinearities in the transfer rate. 
Although these latter two effects could distort the spec- 
trum in a manner qualitatively consistent with experi- 
ment, we find that neither of them plays a significant 
role in these particular experiments. 



aries. These corrugations can be inferred from the un- 
equal spacing of the isodensity contours in Fig. [5J The 
spectrum is a powerful amplifier of these corrugations, 
as they are hardly prominent in the real-space density 



shown in Fig. 4(a) 
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FIG. 4: (a) The density as a function of distance from the 
trap center for a harmonic trap in units of the Thomas-Fermi 
radius £tf, defined as ^tf = y/(fiQ-e )/(m uj 2 /2). The 
density profile corresponds to the LDA contour at t/U a = 
0.018 of Fig. 2(a) The corrugation of the density is observ- 
able, but not dramatic; it would be particularly difficult to im- 
age in experiments looking at columnar integrated densities, 
(b) The hyperfine spectra arising from the density plotted in 
(a). 



III. REFINEMENTS 

As seen previously, GMF+LDA captures the main fea- 
tures of the experimental spectra: sharp peaks occur in 



A. Anharmonicity 

The trap used in the experiments of Ref. Q is a com- 
bination of an Ioffe-Pritchard magnetic trap, which is 
roughly harmonic, and an optical trap, which provide 
highly anisotropic Gaussian confinement. This results in 
a trap with "soft" anharmonicities, increasing the num- 
ber of particles in the low density tails of the cloud. This 
will accentuate the small Su> peaks in the spectrum. The 
presence of anharmonicities is clear in Fig. 4 of Ref. 
where the spatial distribution of the Mott insulator shells 
is far from elliptical. 

We model the trapping potential as 



V tTa . p (x,y,z) 



+ I a (l- e'^l^) +I r (l- e-^+^/^){ll) 

where x lies in the soft "axial" direction while y and z 
constitute the "radial" directions. The 1/e 2 beam waist 
is quoted as 70/im, corresponding to a = 35/im, how- 
ever we find spatial profiles closer to experiment from the 
slightly smaller a = 28/xm and use this value through- 
out. The explicit harmonic terms come from the mag- 
netic trap. The anharmonic Gaussian part has two con- 
tributions, I j from the optical trap and ajVo from the 
optical lattice inducing a further trapping potential, for 
some constants otj, with Ij = I a j + ajVo- The param- 
eters LJ a , u> r , a a , and a r are determined from I 0>a , I o r 
and the quadratic trap frequencies f^vb at Vq = and 
Vq — 4QE Iec by matching the quadratic terms of Eq. pip , 
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FIG. 5: (a) A slice of the anharmonic potential Vtrap(r) 
similar to the one found in experiments, (b) Contour lines 
of constant density in the x-y plane for n = 1, 2, . . . , 5 in the 
anharmonic trap at Vb = 35E ICC . (c) Corresponding spectra 
(using the "spherical trap model" discussed in text) for Vq = 
5, 10, 25, 35£ rcc with TV = 1.4 x 10 7 particles. 
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The quadratic trap frequencies flj.v are given in Rcf. [8] 
as fl r fi — 2tt x 70Hz, fl r _4o = 2tt x 110Hz, Q a ^ — 
2n x 20Hz, and £l a ,40 — 2t\ x 30Hz. The remaining un- 
known parameter I a is chosen to be I _ a = l.l7E Iec so 
that the spatial density profile appears similar to that 
in the experiment. We take I ,r — Io,a though I Q;r 
has little effect on the spatial density profiles. This 
yields u a = 4.8Hz, u r = 67Hz, a a = Q.039E~l, and 
a r = 0.56i?~J to completely characterize the trapping 
potential of Eq. (fTTj) . Note, that while we have chosen 
values to mimic the experimental observations, we have 
not attempted to produce a quantitative "fit" to the ex- 
perimental data. Fig 5(b) shows the isopotential lines of 
our model trap. 

For numerical efficiency, we produce spectra from a 
spherically symmetric model with parameters equal to 
those of the axial direction, which somewhat exaggerates 
the anharmonic effects. As shown in Fig. 5(c) the small 



detuning spectral weight is greatly enhanced by the an- 



harmonicity, producing spectra which are consistent with 
experiments. 



B. Alternative explanations of enhanced 
low-density spectral weight 

Here we examine alternative sources of the enhance- 
ment of the small detuning spectral weight observed in 
experiments. 

Losses. First, we explore the possibility that three- 
body collisions drive the cloud out of equilibrium. Atoms 
are removed preferentially from high density sites, osten- 
sibly enhancing the small-detuning spectral weight. The 
timescale for decay from the n = 5 Mott insulator state 
is 200ms. A characteristic equilibration time is the trap 
period, ~ 10ms. Given the separation of timescales it 
is extremely unlikely that the system is far out of equi- 
librium. Furthermore, the loss rate is effectively zero for 
one- and two-particle site fillings and hence losses are un- 
able to explain the experimentally observed enhancement 
of the n = 1 peak relative to the n = 2 peak. 

Nonlinearities in transfer rate. The probes used 
to measure the spectrum may possibly drive the system 
out of the linear regime where the transfer rate is propor- 
tional to the density. For example, if the transition be- 
comes saturated in the high density regions of the cloud, 
then the observed spectral weight will be reduced at large 
detunings. However, the density dependence of these sat- 
uration effects is slow, making it unlikely that they could 
not be responsible for the dramatic suppression of the ra- 
tio of the spectral weight in the n = 2 and n — 1 peaks. 
A model calculation in the deep Mott regime, where the 
sites decouple, confirms this result. 



IV. SUMMARY 

We calculate the hyperfine spectra of trapped bosonic 
atoms in an optical lattice. We consider the cases of 
harmonic and model anharmonic traps. We show that a 
harmonic trap produces a spectrum which shares quali- 
tative features with the experimental spectra measured 
by Campbell et al. [8j: in the deep superfluid limit one 
has a smooth peak, while in the deep Mott limit, one sees 
several discrete peaks. To reproduce the small-detuning 
spectral weights, however, trap anharmonicities are nec- 
essary. 

We find the spectra are extremely sensitive to density 
corrugation. As an example, the mild density corruga- 
tions which are found in the superfluid near the Mott in- 
sulator boundary are sufficient to produce a pronounced 
splitting of the spectral line. Consequently, the spectra 
are continuous across the superfluid to Mott insulator 
transition. Such continuity is characteristic of a second- 
order phase transition, and makes identifying the super- 
fluid transition difficult. 
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